Spatiotemporal Behavior of Void Collapse in Shocked Solids 
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Molecular dynamics simulations on a three dimensional defective Lennard- Jones solid containing 
a void are performed in order to investigate detailed properties of hot spot generation. In addition 
to the temperature, I monitor the number of energetically colliding particles per unit volume which 
characterizes the intensity of shock-enhanced chemistry. The quantity is found to saturate for 
nanoscale voids and to be maximized after voids have completely collapsed. It makes an apparent 
comparison to the temperature which requires much larger void for the enhancement and becomes 
maximum during the early stage of the collapse. It is also found that the average velocity and the 
temperature of ejected molecules inside a cubic void are enhanced during the collapse because of 
the focusing of momentum and energy towards the center line of a void. 

PACS numbers: 62.50.+p, 47.40.Nm, 82.40.Fp 



In liquid or solid explosives, it is known that hetero- 
geneities such as bubbles or voids decrease a threshold 
shock strength above which a detonation takes place. 
When a shock passes through such defects, they form lo- 
cally excited regions which are reffered to as "hot spots" 
0. Since chemical reactions which turn shocks to det- 
onations are believed to initiate inside such hot spots, 
the mechanism for the generation of hot spots is im- 
portant for the understanding of detonations Q- Since 
shock waves and void collapse are extremely fast micro- 
scopic phenomena, molecular dynamics simulation may 
be a powerful tool for the study of detailed properties of 
hot spot initiation. As well as the studies which treat 
perfect crystals using empirical reactive potentials 
there are a few works on the interaction between shocks 
and defects which can yield hot spots 0, IE IE E] For 
example, Phillips speculated that chemistry would oc- 
cur in a void below a detonation threshold [||. However, 
it is yet to be understood what is essential to the ini- 
tiation of hot spots. Holian et al. recently studied a 
two dimensional Lennard- Jones solid with a planar gap 
which might be regarded as a gap of infinite transverse 
width 8] . They insisted that the essential ingredients to 
hot spot generation are: (i) Enough shock strength for 
the ejection (spallation) of solid molecules into a void, 
(ii) Compression of ejected particles by the far side of a 
void. They estimated the overshooting temperature via 
a straightforward model constructed on the basis of the 
above conjectures. However, since their estimation is not 
quantatively satisfactory (about 3/2 of those obtained 
by the simulations), the validity of their proposal is not 
apparent. Moreover, since void collapse is so fast that 
particles ejected into a void may be far from local equi- 
librium, it is no longer obvious that the "temperature" 
could be a fair index for chemistry. We have to moni- 
tor other quantities to represent the degree of enhanced 
chemistry in this situation. 

In this Letter, I focus on a dynamical aspect of hot 
spot generation to understand what is essential to shock- 
induced chemistry inside a void. Especially, various spa- 



tial profiles of the ejected particles are monitored in order 
to see how the void collapse results in energetic inter- 
molecular collisions which are responsible for chemistry. 
As for the model, I adopt a three dimensional fee crystal 
whose intermolecular potential is the Lennard- Jones one; 



U(r) = 4e 



(7) 



12 



(1) 



whose cut-off length is set to be 4. Oct. Throughout this 
Letter, a and e are set to unity (the LJ unit). Although 
one might wonder about the use of nonreactive potential 
to this kind of issue, this may not be a serious limita- 
tion for the understanding of dynamical aspect of shock- 
induced chemistry; i.e. how the coherent (shock) energy 
relaxes into the thermal energy through the void collapse, 
which is the main aim of this study. The setup of the sys- 
tem is as follows. The axes of x, y, and z are taken to 
coincide with (100), (010), and (001) directions, respec- 
tively. Although there are various methods to realize 
shock waves in molecular dynamics simulations |9j, the 
most popular method is adopted here that a moving pis- 
ton of infinite mass hits the still target. The velocity of 
the piston is denoted by u p . Throughout this study, u p is 
fixed to be 3 which is strong enough to cause spallation. 
Shock propagates along the (001) orientation; i.e. z axis. 
I do not argue orientational dependence which is impor- 
tant to other shock-induced phenomena, such as plastic 
deformation ^(j or phase transitions ^l|. The whole 
system consists of 200 x 200 x 200 unit cells. As usual, 
periodic boundary conditions are applied to x and y di- 
rections. A void consists of blank unit cells, the shape of 
which is cubic or rectangular-solid. Each face is parallel 
to {100}, {010}, and {001}, respectively. Its typical size 
is ranged from 5 to 40 unit cells for each direction. I fix 
the initial temperature to be 0.1 and the corresponding 
density is about 1.053. 

Before proceeding to the results, I make a remark on 
the two kinds of averaging methods adopted in this study. 
One is the average regarding the particles whose initial 
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distance from the upstream wall of the void is less than 
5 (about 3 unit cells). Note that all of these particles 
are eventually ejected into the void after the passage of 
shock front. I have crosschecked the simulations with dif- 
ferent distances from 2 to 15, and the qualitative property 
obtained below is not influenced by a specific choice of 
distance. By this method, we can trace the most ener- 
getic part of the ejected particles. Most of the results 
except for Fig. 21 are based on this averaging method. 
The other method is to take spatial profiles inside a void. 
This method enables us to see where and when the hot 
spot generation takes place. 

Now let us turn to the results of my simulations. First, 
time evolutions of the temperature during void collapse 
are monitored. However, note that the system we are 
considering is noncquilibrium, where the definitions of 
thermodynamic quantities including the temperature are 
precarious [l2|. Here 1 adopt the variance of z velocity 
as the "temperature" , which is denoted by T z . For com- 
parison, I also monitor the variance of x velocity which 
is denoted by T x . Note that the variances can be dif- 
ferent values for different directions; e.g. T x is smaller 
than T z in this situation. Temporal behaviors of T z , T Xl 
and the average z velocity denoted by v z are shown in 
Fig'H Spallation occurs at t = 5 where v z suddenly in- 
creases up to about 2u p = 6 (the spall velocity). Af- 
ter the spallation, all of the quantities increase gradually 
until the ejected molecules collide with the far side of 
the void (approximately t — 10). It should be stressed 
that v z increases up to about 2.6u p , which means that 
there is some kind of momentum transfer into the ejecta 
gas. (We see this acceleration by observing the velocity 
distribution functions in the next paragraph.) Also it 
should be remarked that the maxima of T z and v z are re- 
alized simulataneously, which is just before the collision 
of ejecta gas and the downstream wall of the void. It is 
not only the compression by the void wall that is respon- 
sible for heating of the ejecta gas but also some kind of 
focusing of momentum and energy. However, recall that 
the above results only concern with the sample molecules 
which constitute the proceeding part of ejecta gas. As for 
the whole gas, observation of the spatial distributions in- 
side the void is needed, which we will see later in detail. 

To see the above circumstances regarding the accelera- 
tion of v z in detail, velocity distribution function is shown 
in Fig. [3] Strong displacement from the Maxwellian are 
seen after the passage of the shock front. Note that a 
distribution function of v x is also observed to show an 
exponential distribution just after the spallation, which 
then eventually relaxes to the Maxwellian. On the other 
hand, we can see that the distribution function of v z has 
a widely stretched tail for the positive direction, which 
develops until the collision with the far side of the void. 
Note that the tail by far exceeds the spall velocity. It is 
obviously responsible for the increase of v z , since the tail 
develops only towards the positive direction. The devel- 
oping tail means that there is a momentum transfer into 




FIG. 1: Temporal behaviors of the longitudinal (solid line) 
and the transverse (dotted line) temperatures, the average z 
velocity (dashed line). The number of energetically colliding 
particles which is introduced afterwards is also shown by the 
dot-dashed line (scaled by 10). The sample molecules are 
those initially locate near the upstream wall of the void. (See 
text.) The void initially consists of 30 x 30 x 30 blank unit 
cells. 
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FIG. 2: Velocity distributions of v z . Time ordering is: (i) solid 
lines (t = 0; before the spallation), (ii) dashed lines (t — 6.0; 
just after the spallation), and (iii) dotted lines (t = 9.5; just 
before the collision with the far side of the void). Note that v z 
is shifted by — u p . All of the parameters including the sample 
molecules are the same as in Fig. 



the ejecta gas, which may be attributed to the focusing. 
As is visualized in Fig. [3| since the shock pressure also 
pushes the side wall of the void to form a parabolic sur- 
face, the momenta of the collapsing wall can focus into 
the center line of a void, which eventually causes jetting 
and accelerates v z . Note that the acceleration during the 
flight does not happen in a rectangular-solid void with 
sufficiently large longitudinal length where the parabolic 
surface is not formed. It is obvious that larger v z yields 
more energetic collision, which eventually leads to higher 
overshooting temperature inside a hot spot. Phillips et 
al. |5j also observed deviations from equilibrium distri- 
butions in a collapsing void using a two dimensional LJ 
crystal. However, their observation focused only on a 
relaxation to equilibrium, whereas the transition from 
equilibrium to noncquilibrium is traced here. Also, the 
positive tail which prevails spall velocity are not seen 
there, which is the original contribution of this work. 
Apart from the temperature and the spall velocity, we 
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FIG. 3: (a) A slice-cut view (parallel to {100} plane) of the 
collapsing void and (b) its velocity vector field, (t — 9.8.) A 
shock propagates upwards. 

then turn to the enhancement of chemical reactions which 
is our original concern regarding hot spots. Since ejected 
particles are strongly nonequilibrium as we have seen 
above, it is no longer expected that the "temperature" 
T z can be a fair index of chemistry. Although chemistry 
involves many intramolecular degrees of freedom, strong 
intermolecular collision is essential to the excitation of 
those degrees of freedom. Therefore it is quite natural 
and direct to count intermolecular collisions for the es- 
timation of chemistry. In this study, I count the cou- 
ple of particles whose distance is smaller than a certain 
length, r c . Then the number of the colliding particles 
is normalized by an initial volume of a void. Hereafter 
this intensive quantity is denoted by N c . The value of 
r c can be various according to the molecules we are to 
consider; in this case, I choose r c = 0.88 which corre- 
sponds to U(r c ) — 2.5. Note that U(r c ) gives an energy 
threshold for a reaction to occur. However, since LJ po- 
tential is adopted here, we do not have the quantitative 
validity for the choice of r c . Therefore, I crosscheck the 
behaviors of N c for several values of r c and no qualita- 
tive difference is seen regarding the results obtained in 
this Letter. Typical temporal behavior of N c is shown 
in Fig. ^ Note that N c is zero until the void collapse, 
which means that strong intermolecular collisions take 
place only in the void. Also it should be remarked that 
the peaks of N c and of T z are not simultaneous; It takes 
a certain amount of time that N c becomes maximum af- 
ter the peak of T z . Indeed, this delay means that en- 
hanced chemistry needs higher density than overheating; 
Namely, the delay is a relaxation time for the density to 
increase in the collapsing void. To confirm this idea, it is 
observed that the delay becomes longer for larger voids. 
This situation becomes clear by monitoring spatiotem- 
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FIG. 4: Spatial profiles of (a) T z and (b) N c represented by 
thick lines together with the density profiles (dotted lines) 
which are rescaled suitably. (Multiplied by (a) 20 and (b) 
0.25.) Each graph is a snapshot at the moment of the maxi- 
mization of each quantity: (a) t = 10.0 for T z and (b) t = 11.0 
for N c . The initial density profiles are also imposed (thin 
lines). 



poral behavior of the density profile inside the void. Fig 
01 shows the spatial profiles of N c , T z , and the density, 
through which we can see in datail where and how the 
energetic molecular collisions and the overheating take 
place. Note that these two snapshots are taken at the 
moments when T z and N c become maxima, respectively. 
The maximum N c is realized after the whole region of 
the void has collapsed, while the preceeding maximum 
temperature takes place in the relatively dilute region on 
the far side of the collapsing void. It is important to no- 
tice that N c becomes maximum after the void has totally 
collapsed, whereas T z is maximized at the beginning of 
collapse. The total collapse yields high density to en- 
hance energetic intermolecular collisions. Therefore the 
time lag between the maxima of T z and N c are the time 
needed for the total collapse of the void. 

The discussions so far finish spatiotemporal behaviors 
of the ejecta gas and enhanced chemistry. As the final 
part of my results, I show how the size and the shape of 
voids affect chemistry. Although voids in real materials 
can take various forms, it may be decomposed into two 
essential ingredients: the longitudinal length and the nor- 
mal cross section. In order to specify which component is 
important to the initiation of chemical reactions, I check 
rectangular-solid voids changing the transverse and the 
longitudinal lengths independently. First I show the lon- 
gitudinal length dependence with the transverse length 
(and hence the cross section) fixed. We can see that N c 
decreases for larger longitudinal length, whereas T z in- 
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FIG. 5: (a) The longitudinal and (b) the transverse lengths 
dependences of T z (filled circles) and N c (blank squares, mul- 
tiplied by 20). The length in the fixed direction of the void is 
30 in each case. 

creases. As for T z , this result is consistent with Q where 
a considerable longitudinal length was needed for the en- 
hancement of temperature. However, the decrease of N c 
again shows that the temperature cannot be a reliable 
index for enhanced chemistry in this nonequilibrium sit- 
uation. On the other hand, the cross section dependence 
of T z and N c seems quite different from the longitudinal 
length dependence. The temperature T z is maximized 
when the void is cubic, whereas N c increases to saturate 
for the larger cross section. It should be stressed that 
the cross section of 45 x 45 (approximately 30 x 30 unit 
cells) is enough for the saturation. The maximization 



of T z with a cubic void may be a result of focusing dis- 
cussed in the context of the acceleration of v z . Also note 
that a larger cubic void can yield higher temperature and 
v z . The saturation of N c for larger cross section means 
that the focusing in a cubic void does not seriously in- 
fluence the number of energetic intermolecular collisions. 
Recalling that higher density (the total collapse) is also 
required for the increase of N c , it is plausible that the 
intensity of void collapse is not affected by the focusing, 
whereas the foregoing part of ejecta gas is accelerated 
to enhance the temperature. Hence it is concluded that 
N c is rather insensitive to the shapes of voids, whereas 
T z has the apparent tendencies. In any case, nanoscale 
voids (approximately 30 unit cells in each dimension) is 
enough for the enhancement of energetic intermolecular 
collisions. 

To summarize, through the observations of spatiotem- 
poral evolutions of temperature and density inside a void, 
we have seen that energetic intermolecular collisions are 
most likely to occur when a void is totally collapsed. 
This is not simultaneous with the maximization of tem- 
perature which is realized at the beginning of the col- 
lapse. This time lag is due to the difference between 
the optimum densities where the each quantity is maxi- 
mized. The shape of void affects the overshooting tem- 
perature, but cannot seriously affect the number of en- 
ergetically colliding particles per volume. Especially, the 
intensity of enhanced chemistry is found to saturate for 
nanoscale void whereas the enhancement of temperature 
needs much larger void. Since chemical reactions origi- 
nate from intermolecular collisions, monitoring only the 
temperature is not enough to explore the various proper- 
ties of hot spot generation, such as the void-shape depen- 
dence on the enhanced chemistry. Also, velocity distribu- 
tion functions are found to deviate from the Maxwellian. 
These strong nonequilibrium effects must be taken into 
account for the analysis of spallation phenomena. 

The author is grateful to H. Kaburaki and F. Shimizu 
for helpful discussions and encouragements. 
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